Construction of exact solutions of Bloch-Maxwell equation based 

on Darboux transformation 

Maciej KunaQ 
Abstract 

A new strategy, using Darboux transformations, of finding self-switching solutions of ip = 
[H,f(p)] is introduced. Unlike the previous ones, working for any / but for Hamiltonians whose 
spectrum contains at least three equally spaced eigenvalues, the strategy does not impose any 
restriction on the discrete part of spectrum of H. The strategy is applied to the Bloch-Maxwell 
system. 
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I. INTRODUCTION 

The nonlinear von Neumann equation (NvNE) 

V =[#./(/>)] (1) 
was recently shown to be relevant for two apparently disconnected fields of research: Optical 
solitons 5] and chemical kinetics f|. In both cases we are interested in solutions p that have 
a form of n x n matrices. In optical applications the dimension n is related to n-level atoms. 
For chemical-type kinetics the kinetic variables are related to real and imaginary parts of off- 
diagonal matrix elements of p evaluated in the basis of eigenvectors of H. Since any chemical 
kinetics involves a finite number of substrates, the dimension n must also be finite. It has 
to be stressed, however, that a class of infinite-dimensional solutions of JI} was also found, 
albeit for a rather artificial example, in Moreover, as shown in []| the NvNE contains as 
special cases many known lattice equations that also lead to infinitely-dimensional solutions. 

The soliton-type solutions of (JTJ were first found in |2j for f(p) = p 2 and a special class 
of Hamiltonians H . The Hamiltonians H could be characterized as those that contained 
at least three equally spaced eigenvalues in discrete parts of their spectra. There are, of 
course, many physical examples of such systems, the most notable cases including harmonic 
oscillators, quantum fields or spin systems. However, the problem that was not solved until 
now was how to construct solutions for H that has an arbitrary spectrum. In the present 
paper we show how to deal with any discrete spectrum and illustrate the technique on a 
harmonic oscillator (HO) and a hydrogen- type atom (HA). 

The method we employ is based on a Darboux transformation. The main advantage of 
the technique is that it allows to construct solutions that are practically impossible to find 
by other means, especially if n is large or infinite. The disadvantage lies in the need of 
starting with a known solution that, in addition, guarantees nontriviality of the solution 
generated by the Darboux transformation. 

We work out explicitly two important cases: n = 2 and n — 3. The cases are important 
for two main reasons. Firstly, the two dimensional case is easy to handle and allows for a 
seed solution having a sufficient number of arbitrary parameters, no matter what spectrum 
of H one encounters. Secondly, once we have two-dimensional seed solutions, we can use 
them to construct highly nontrivial "self-scattering" solutions involving any even n. If n is 
odd we can employ a combination of two- and three-dimensional seed solutions. 



In the present paper we restrict the analysis to the essential building blocks with n = 2 and 
n = 3. But even these cases, especially n = 3, are interesting in themselves. For quantum 
optics we can apply them to three-level atoms and the result is a dynamics analogous to 
the celebrated "sech" soliton In application to chemical kinetics the case results in an 
Oregonator-type dynamics with four substrates 

In both cases it is sufficient to work with quadratic nonlinearities. Kinetic autocatalytic- 
type equations are found immediately by considering separately real and imaginary parts of 
matrix elements (k\p\l) where \k), \l) are eigenvectors of H normalized to (k\l) = 5m- 

In order to relate the dynamics to n-level atoms one makes the trick introduced in j^j]. 
Take a solution p of ip = [H, p 2 } and employ the algebraic formula [H, p 2 } = [Hp + pH,p}. 
For a given solution p rewrite Hp + pH as H' + V where H' has the same form as H but 
with modified parameters. Then V is a time-dependent interaction term that has a form 
—d ■ E(t, 0) where E(t, 0) is a field typical of an optical soliton at the origin x = 0. As 
pointed out in [^| the construction is analogous to SUSY quantum mechanics since we have 
in addition the Darboux transformation p f— > p\ that allows for the transformation V i— > V\. 
We leave these details for a later work but in the present paper concentrate on technicalities 
leading to the solutions for arbitrary H . 

II. LAX PAIR AND DARBOUX TRANSFORMATION 

Following one starts with three overdetermined linear systems (Lax pairs) 



where p, H are self-adjoint operators on some Hilbert space 7i, pk> z k are complex numbers, 
/ is real function and \iff.) are elements of 7i for k = 1, 2, 3. If we fix an operator H to be 
time-independent then there exist an operator p and a vector | that fulfil the Lax pair 
only if ip\ifik) = [H, f(p)]\(fk}- The basic theorem jj| is the following: 

Theorem. Assume \ipk) for k — 1,2,3 are solutions of (0 EJ) respectively and |^[1]), p[l], 
are defined by 



ZkWk) = (p ~ p k H)\ip k ), 

#*) = — f{p)Wk)i 
Pk 



(2) 
(3) 




(4) 
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Then 



rfl] = ( 1 + «Z»p) <I ( 1 + »Z»p) 1 (5) 



/(P[l]) = (1 + ^^P)/(P)(1 + ^^P), (6) 

v /i 2 ' v Pz ' 

p = MM. ( 7) 

mm) 



Zl \m) = ( P [i]*-» k H)\m), (8) 

mi) = -f(p[MT\m), (9) 
pi 

If we put p 3 = p 2 and | v^i ) = Iv^) then p[l] is also a density matrix and spectra of p 
and p[l] are identical. In this case © is simply equivalent with definition of function / of 
operator p[l] given by the spectral theorem. We can proof jll| that if p is a solution of (JTJ) 
then p[l] is also a solution of (JTJ). 



III. TWO-DIMENSIONAL CASE 

We start with a Hamiltonian having discrete spectrum: H = J2n h n Qn, where Q n and 
h n are spectral projectors and eigenvalues of H respectively. Next we take a density matrix 
p with two spectral projectors Pi and P 2 which satisfy the following conditions: [Pi,Q n ] = 
[P2,Qn] = for all n ^ 1,2 Therefore P = Pi + P 2 = Qi + Q2 is a projector on the 
two-dimensional common eigensubspace. Let us define p\ = PpP, Pi = PPP, p' — (1 — 
P)p(l— P) and H' = (1— P)P(1— P). Then equation ip = [H, f(p)} separates into two pieces 
ip\ = [Pi,/(pi)] and ip' = [H',f(p')}. By using the relation [P, Pi] = and properties of 
such evolution P = the first equation can be described by spectral projectors: 

ip\ = z(A x Pi + A 2 P 2 ) = i(Ai - A 2 )A 
[HiJ(pi)] = [(/(Ai)-/(A 2 ))P 1 ,P 1 ], (10) 

so 

iP 1 = [a(X)H 1 ,P 1 ) 

and similarly: 

iP 2 =[a(X)H 1 ,P 2 ], 
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where a (A) = ^ A ^_{|, A2 ' > ■ Solutions of these equations are the following: 

Pi(t) = e -^WHitp^ e ia(X)H lt 

for i — 1, 2 and we have the general solution of two-dimensional problem: 

Pl (t) = e- ia{x)Hlt Pl (0)e iaWHlt . 

We want to use it in further constructions, therefore we have to examine its properties. 
Let us assume that p\ together with H\ and fixed z^, p and \ipi) fulfil the Lax pair on PTi: 



= (Pi ~ ^ H i)\fi) 

1 



i\<Pi) = -fc(f>i)\<Pi), (11) 



where / c (pi) = /(pi) + cP. 
Put = e ia(A)Hl *|(^i), then 



z^i) = Oi(o)-/^iM> 

= {-/ c (pi(0))-a(A)F 1 }|^ 1 ). (12) 

So IV'i(O) — e~ t ^^ pl ^~ a ^ Hl ^ t \'ipi(0)). It is easy to check that the generator of time 
evolution of commutes with pi(0) — pB.\. 
Taking this into account we obtain 

^|^i(0)) = ( Pl (0)- fiHJIMO)) 

\ Vl (t)) = e -^W^* e - i ^/°(«(°))- a ( A )^}*|^ 1 (0)). (13) 

We can conclude that if p\ is a solution of the two-dimensional problem then it is enough 
to consider only the Lax pair for time t = 0. 



IV. THREE-DIMENSIONAL CASE 



We construct a density matrix p using only three eigenprojectors: P±, P2 and third 
spectral projector is equal to Q3. So p' = A3P3 = X3Q3. Denote by |cuj(0)) eigenvectors of 
p(0) , so Pi(0) = |u;j(0))(u;j(0)| and denote by Pij(0) = \u)i(0)){(jjj(0)\. Now we want to add 
to the Lax pair (|TT|) a component on this one-dimensional space in such a way that together 
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they will be the three-dimensional Lax pair with slightly modified /. For this case we can 
also restrict to initial conditions and for given z^, fi and I13 we can choose A3 satisfying the 
Lax pair: 

z> 3 (0)) = (A 3 -M 3 )k 3 (0)) 

\u 3 (t)) = M0)>. (14) 

Put \<f(t)) = 7i|<^i(i)) + 73)^3(0)) for arbitrary 71,73 G C. Then we have the three- 
dimensional Lax pair: 

= ip-V H )\v) 

iW) = kf( P )-f(\s)M = \f{p)\<p) (15) 
and we have to consider only 

z^(0)) = (p(0)-A*H)|^(0)> 

\(p(t)) = e - ia ^ m e- i{ ^ f ' ipm - aWH}t \ij(0)) = U\ij(0)) (16) 

with \ip) = e ia ^ Ht \v). 

Now we have the Lax pair expressed only by initial conditions and can construct the 
operator P which defines the Darboux transformation (see[7j): 

\<p(t))(ip(t)\ _ U\^(0))(^(0)\e i{ ^ f ' {pm - a{x)H}t e ia ^ m 

p " (<p(t)Mt)) ~ ¥(MpW) 

= e -ia(X)Ht p . nte iaWHt = Vp . ntV * y ( 17 ) 



where 



Then 



p \mm\ _ \mm\ 
mt (my®) ~ (m\m)' 



P [i](t) = (i + ^P) P (t)(i+^P) 

= ^(1 + ^P mi )p(0)(l + ^—^PintW* = Vp mt [1] (t)V*. (18) 
From the first equation of Lax pair (J2J) we have 

H\^} = -(p(0) - z,M) 
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so 



( V(p(0)) - a(X)HM = V( P (0)) - a(A)KO) - ^])|^> 



and the second equation of Lax pair reads 



#> = V(p(0))-a(A)[p(0)-^])|^. 



It can be written as 



|^ = e f{/'«0))-«(AM0)-z M ]}t|^^^ 

= e f{Et 1 ^m}*| V , (0))) (19) 
where = f'(\i) — a(A)[Aj — z M ]. Inserting this into the definition of Pj nt we obtain 



p -s^ x r ( "^ (o,i#)iwoiM 



where a 4J = (^(0)|^(0))(^(0)|^(0)), b tj = - £&2. Since 6 n = b 22 = b 12 = b 21 = b and 
(3(X 3 ) = -a(\)nh 3 then 633 = 0. So F(t) = £? =1 a u e~ ibiit = e- lbt (a n + a 22 ) + a 33 - If 73 = 
then the projector P int does not depend on time, so we need at least a three dimensional 
space to have a nontrivial Darboux transformation. We are in position to give the form of 
Pint [1] (t) in terms of spectral projectors of p(0): 

Pmt [l]{t) = (l + ^ Pint )p(0)(l + ^^P int ) 

= fW) + fe4{E^(«) + E%^(o)}, (21) 

^ l r J IpI i=i 

where q = (/x — /l)ajje~* 6 "*[(Ai — A,k/,,< ''' ' + (Ai — A fc )a fcfc e~ l6fcfc *] and c^- = aije~ lbijt [(\j — 
\i)(na ii e~ lb " t +~pajje~' lbj: > t ) + ((~p — pi)\ k + —~p\i)a kk e~ lbkkt ] and here j, in the first term 
and k in the second term mean the remaining indexes. 

V. MODELS 

We want to see whether we can find appropriate /i and a density matrix p for a construc- 
tion of "self-scattering" solution for any Hamiltonian H with a discrete part of spectrum 
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and nonlinearity /. We use the spectral representation of H and take any three of its 
eigenvectors: 



' h x o o\ 
H= h 2 
v o h 3 J 

We can decompose H in two pieces in two different ways: 

/ Uh - h 2 ) o o N 



(22) 



H = 



V 



o 
o 



-\{h 1 -h 2 ) 







h^-\{h 1 + h 2 )j 



+ [h l --(h 1 -h 2 )]I 



(23) 



for /i 3 — \{hi + ^2) 7^ and 
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\ 



V 



+ a/ 



(24) 



35 
2B 

for equally-spaced spectrum. Clearly, if the density matrix p is a solution of the equation 
ip = [H,f(p)], then it also satisfies the equation with H' = H + CI, so we can consider 
Hamiltonians with special forms. In the first case, for the operator 

( Qi-^{hi-h 2 ) q N 



p- pH 



V 



Q 




Q2 + 1(hi - h 2 ) 







(25) 



g 3 - f/[/i 3 - + h 2 )\ ) 

we can choose Q\ — Q 2 = a(hi — h 2 ), (3 2 (hi — h 2 ) 2 )A\g\ 2 ) \g\ 2 = (3 2 [h 3 (hi + h 2 — /i 3 ) — hih 2 ] 
and 0i + g 2 = 2g 3 — 2ah 3 + a{h\ + h 2 ) , where p — a + i(3 to have an appropriate eigenvector. 
The compatibility of the third assumption is satisfied if we arrange eigenvalues of H in 
such a way that h 3 is a number between h\ and /12. Under these assumptions we can take 
g l = a (h l - h 3 ) + 03, 02 = cn(/i 2 _ ^3) + Q3 and I0I 2 = /3 2 [h 3 (hi + h 2 - h 3 ) - h x h 2 \. To get 
p positive we need 0«)O and 0102) |f?| 2 - We can easy control these conditions by choosing an 
appropriate 3 , a and (3. In order to have p as a density matrix we have to put a = hl +^% h3 ■ 
We will latter discuss this condition. For such a specification p — pH has the eigenvalue 



z. 



03 + |(1 - 30 3 ) - i(3[h 3 - i(/n + h 2 )\. 



M ^3 1 2 

In the second case: 



p'-pH 



q'i - liB 0' 



Q' 






3^5 
3 - 2^5 y 



(26) 
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has an appropriate two-dimensional eigenspace if, for example, we take g[ — g' 2 = —2aB 
and \g'\ 2 = j3 2 B 2 . To obtain density matrix we can put g' 3 = |, g[ = | — aB, g' 2 = | + «B, 



>'| 2 = /3 2 5 2 and l^l 2 ^ 2 ^. Then z ll = \- 2/iB. 
Let us choose a hydrogen-like atom Hamitonian (HA) as an example of an inhomogeneous 
spectrum and a harmonic oscillator (HO) for the homogeneous one. We choose neighbouring 
levels for simplicity: h[ HA) = ~ B ^ A) , h^ 1 ^ = ~£+ 2 )* ' ^3*^ = "tJttJ^ • ^ nave som tio n s 
not so complicated in form and more legible we take a — 0. Otherwise the calculations 
become very complicated. We find g ( 3 HA) = \, g[ HA) = (n+2) ^"^^fllf" +24n+8) , g ( 2 HA) = 



i (6ra 3 +15n 2 +12n+4)-2n 3 -3n 2 



and |f? 



6n 2 +12ri+4 

obtain the same type of solution 



(HA) 1 2 



/3 2 (g( gA )) 2 (4n 2 +8n+3) 
n 2 (n+2) 2 (rt+l) 4 



6n 2 +12n+4 

. Then for both Hamiltonians we 



Pint[L](t) 



1/3 + C(*) 








\ 



v m 



1/3 - c(o 
m 1/3 / 



(27) 



For HO we have 



C(t) = \g'\(i 



2\ e 



'l 2 t ) 



2+ LD \ 2 e~ 



\gft 



D\ 

35' 



= \g'\ tanh(^- + ln(^)) = \g'\ tanh(0'i + 0') 



2\g'\D(l + i)e~ 



2 + \D\ 2 e 
D 



2|o'| 2 t 



\g'\H 



V2\D\ 



\D 
V2 



^'|sech(^- + ln(^i)) = Z'sech(6't + 



(28) 



and for HA 



2d(n) + \D\ 2 e 3 

= | g ( gA >ltaDh( |g(H ^ )|2 * +ln( -j£L)) = |^)|tanh(^t + ^>) 

i e (^)| 2 t 

m = -— — ^'"'"i'lL, =- r::z v * A vch(^r-^+M^==)) 



2\g {HA ^\e(n)De — 3 — De(n) . (HM . , Ag {HA) \H |D| 

— 2 , (HA) = V ; \g {HA) \sech C 1 + ln(— 1— I 

2d(w) + [£>l 2 e ^/2d(n)|D| P \J^d{n) 



= Z^sech^ HA H + ^), (29) 

where d(n) = 4(2n + l)(n+ l) 3 , e(n) = (2n + 1) (n + 2) + in^/(2n + 3) (2n + 1) and D = ffi. 
Rewriting pj n t [!](£) in spectral basis of the Hamiltonian and renumerating in the order of 
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increasing eigenvalues we find: 



( 1/3 y/mm) m N 

Pint[l}(t)= V2ft(Z(t)) 1/3 -*V2»(^(*)) • (30) 
^ C(t) iv^^t)) 1/3 y 

If we compare this expression with formula (16) in [5] we can see that we have two differences. 
First, the diagonal elements are different but in both cases constant. The second difference is 
more important. (16) describes p[l](i) but not ^ nt [l](t) and we have to sandwich the latter 
between e ia \ x ) Ht and its inverse in sense of f|18|) to obtain the former. After this operation 
the symmetry of p[l](£) for HO and HA will be slightly different. For HO we have exactly 
the same symmetry as in p| and for HA £ in the first row is multiplied by e ia{X){hl - h2)t and 
in second row by e ia ( x )( h 2-h 3 )t ^ QT j^q are ^he same ). 

Obviously these differences do not disturb the idea of p| to construct a time dependent 
Hamiltonian describing a three-level perturbation of not only HO but for any system, e.g. 
HA interacting with an optical soliton. The only modification comes from non-equal spacings 
between levels of the system and produces appropriate frequencies of oscillation. We can 
easily observe that the only modification of H is in multiplication by the factor | that, for 
the quadratic case, is equal exactly to a (A). Therefore we can start with Hamiltonian times 
| and get an appropriate H a . 

The relation ip = [H, p 2 ] = [Hp + pH, p] = [h, p] is given up to two parameters E\ and e 2 : 
h = (H + £i)p + (H + Ex)p + e a l. Putting e Y = -Ma an d £2 



2^21 we get: 



/ 2/3/ii \ 

2/3h 2 
\ 2/3h 3 J 



( a(t) ^ 
o(t) b(t) 



H n + V = H n -d-E, 



(31) 



V b(t) J 



where a{t) = (h 2 - /i 3 )v / 23?(£(t))e i2/3(fel - ?i2) * and b(t) = -i(h 2 - h)^^))^ 2 ^ 3 '^ ■ 
From this we can see that: 



E = (E x , E y , E z ) = (E x (0)sech(9t + i9) cos(cjt), E y (0)sech(6t + •&) sin(wt), 0) (32) 

and 
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d, 



( 







d(t) cos(c^) \ 
a(t) cos(ui) b(t) cos(cot) 




b(t) cos(ujt) 







, (iy 







d(t) sin(wt) \ 
a(i) sin(u;£) sin(a;t) 
6(t)sin(u;f) j 



(33) 

a(t) = (h 2 - h 3 )^/2^{Z)e l2 / 3 ^~ h ^, bit) = -i(h 2 - h 1 )y/2$f(Z)e i2 M ha - h *) t . Hence we get 
potential V in shape of well-known McCall-Hahn "sech" soliton jjj. Another choice of 
parameters £j leads us to a very wide class of electromagnetic impulses We can use 

now this solution to construct another one for the same potential in accordance with the 
SUSY scheme. The possibility follows from the richness of the set of solutions we obtain by 
the Darboux transformation. 

Now we are in position to check that we have not only the solution of the Bloch equation 
but also of the Maxwell one. First we calculate the components of atomic polarization: 



P x = Tr(p4) = 4sech(0t + #) cos(cut)[(h 2 - h 3 )^(Z) 2 + (h 2 - h^^Z) 2 } 

P y = Tr(pdy) = 4sech(0t + 0) sin(cjt) [(h 2 - h 3 )$l(Z) 2 + (h 2 - h^Z) 2 } (34) 

Let us remember that (h 2 — h 3 ) and (h 2 — hi) have the same sign. 

Because McCall-Hahn is a steady-state pulse we can use coordinate <j = t — -, where v is 
the constant pulse velocity. Hence we have i— > ^ and i— > — ^ then Maxwell equation 
has the following shape: 

d 2 - 47T d 2 - 

By comparing ()32j) and ()34j) we can see that Maxwell equation is fulfilled if we choose 
parameters in such a way: 

16n[(h 2 - h 3 MZ) 2 + (h 2 - h)^Z) 2 ] = ((^) 2 - l)J5 o (0), 

where o means x or y. 



VI. COMMENTS 



We have described a construction of solutions for the two main "building blocks" : 2 x 2 
and 3x3. How to get from 2x2 blocks to infinite-dimensional and irreducible solutions 
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was discussed in |4j. However, in infinite dimensional examples the issues of normalization 
to Tr p = 1 involve certain subtleties that require a separate treatment. These difficulties 
do not occur if one restricts the discussion to arbitrarily large but finite matrix dimensions. 
The method we have described in the present paper allows for immediate generalizations to 
arbitrary finite dimensions if one follows the strategy employed in ^j]. 

The technique gives a method of constructing three-dimensional solutions appropriate for 
different situations. It has to be stressed that not only can we use any Hamiltonian and any 
type of nonlinearity, but it seems that we can get arbitrary values for the eigenvalues of p for 
a homogeneous nonlinearity. To see this it is enough to observe that we can generalize the 
scaling and the shifting properties which were observed in [2j the for quadratic nonlinearity. 
If g is solution of NvNE for f(x) = x k then fy^( ( Tr p)fc-i ) also satisfies it and p + a when a = 
is solution if we put g(x) = f(x + a) instead of f(x). Because Qi and g 2 are very simple 
functions of the parameter a and g 3 , and g depends only on /3, we can produce positive 
"predensity" matrices with wide spectrum of eigenvectors. How "big" is this set of solutions 
in the set of all the solutions of NvNE will be considered elsewhere. 
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